Investigating the potential of single-cell DNA methylation data to detect allele-specific methylation and imprinting.

Allele-specific methylation (ASM) is an epigenetic modification whereby one parental allele becomes methylated and the other unmethylated at a specific locus. ASM is most often driven by the presence of nearby heterozygous variants that influence methylation, but also occurs somatically in the context of genomic imprinting. In this study, we investigate ASM using publicly available single-cell reduced representation bisulfite sequencing (scRRBS) data on 608 B cells sampled from six healthy B cell samples and 1,230 cells from 11 chronic lymphocytic leukemia (CLL) samples. We developed a likelihood-based criterion to test whether a CpG exhibited ASM, based on the distributions of methylated and unmethylated reads both within and across cells. Applying our likelihood ratio test, 65,998 CpG sites exhibited ASM in healthy B cell samples according to a Bonferroni criterion (p < 8.4 × 10-9), and 32,862 CpG sites exhibited ASM in CLL samples (p < 8.5 × 10-9). We also called ASM at the sample level. To evaluate the accuracy of our method, we called heterozygous variants from the scRRBS data, which enabled variant-based calls of ASM within each cell. Comparing sample-level ASM calls to the variant-based measures of ASM, we observed a positive predictive value of 76%-100% across samples. We observed high concordance of ASM across samples and an overrepresentation of ASM in previously reported imprinted genes and genes with imprinting binding motifs. Our study demonstrates that single-cell bisulfite sequencing is a potentially powerful tool to investigate ASM, especially as studies expand to increase the number of samples and cells sequenced.

A comparative atlas of single-cell chromatin accessibility in the human brain.

Recent advances in single-cell transcriptomics have illuminated the diverse neuronal and glial cell types within the human brain. However, the regulatory programs governing cell identity and function remain unclear. Using a single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), we explored open chromatin landscapes across 1.1 million cells in 42 brain regions from three adults. Integrating this data unveiled 107 distinct cell types and their specific utilization of 544,735 candidate cis-regulatory DNA elements (cCREs) in the human genome. Nearly a third of the cCREs demonstrated conservation and chromatin accessibility in the mouse brain cells. We reveal strong links between specific brain cell types and neuropsychiatric disorders including schizophrenia, bipolar disorder, Alzheimer's disease (AD), and major depression, and have developed deep learning models to predict the regulatory roles of noncoding risk variants in these disorders.

Single-cell DNA methylation and 3D genome architecture in the human brain.

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Differential DNA methylation and changing cell-type proportions as fibrotic stage progresses in NAFLD.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is characterized by changes in cell composition that occur throughout disease pathogenesis, which includes the development of fibrosis in a subset of patients. DNA methylation (DNAm) is a plausible mechanism underlying these shifts, considering that DNAm profiles differ across tissues and cell types, and DNAm may play a role in cell-type differentiation. Previous work investigating the relationship between DNAm and fibrosis in NAFLD has been limited by sample size and the number of CpG sites interrogated. RESULTS: Here, we performed an epigenome-wide analysis using Infinium MethylationEPIC array data from 325 individuals with NAFLD, including 119 with severe fibrosis and 206 with no histological evidence of fibrosis. After adjustment for latent confounders, we identified 7 CpG sites whose DNAm associated with fibrosis (p < 5.96 × 10-8). Analysis of RNA-seq data collected from a subset of individuals (N = 56) revealed that gene expression at 288 genes associated with DNAm at one or more of the 7 fibrosis-related CpGs. DNAm-based estimates of cell-type proportions showed that estimated proportions of natural killer cells increased, while epithelial cell proportions decreased with disease stage. Finally, we used an elastic net regression model to assess DNAm as a biomarker of fibrotic stage and found that our model predicted fibrosis with a sensitivity of 0.93 and provided information beyond a model based solely on cell-type proportions. CONCLUSION: These findings are consistent with DNAm as a mechanism underpinning or marking fibrosis-related shifts in cell composition and demonstrate the potential of DNAm as a possible biomarker of NAFLD fibrosis.

Ten-eleven translocation protein 1 modulates medulloblastoma progression.

BACKGROUND: Medulloblastoma (MB) is the most common malignant pediatric brain tumor that originates in the cerebellum and brainstem. Frequent somatic mutations and deregulated expression of epigenetic regulators in MB highlight the substantial role of epigenetic alterations. 5-hydroxymethylcytosine (5hmC) is a highly abundant cytosine modification in the developing cerebellum and is regulated by ten-eleven translocation (TET) enzymes. RESULTS: We investigate the alterations of 5hmC and TET enzymes in MB and their significance to cerebellar cancer formation. We show total abundance of 5hmC is reduced in MB, but identify significant enrichment of MB-specific 5hmC marks at regulatory regions of genes implicated in stem-like properties and Nanog-binding motifs. While TET1 and TET2 levels are high in MBs, only knockout of Tet1 in the smoothened (SmoA1) mouse model attenuates uncontrolled proliferation, leading to a favorable prognosis. The pharmacological Tet1 inhibition reduces cell viability and platelet-derived growth factor signaling pathway-associated genes. CONCLUSIONS: These results together suggest a potential key role of 5hmC and indicate an oncogenic nature for TET1 in MB tumorigenesis, suggesting it as a potential therapeutic target for MBs.

Age-related DNA hydroxymethylation is enriched for gene expression and immune system processes in human peripheral blood.

DNA methylation (DNAm) has a well-established association with age in many tissues, including peripheral blood mononuclear cells (PBMCs). Compared to DNAm, the closely related epigenetic modification known as DNA hydroxymethylation (DNAhm) was much more recently discovered in mammals. Preliminary investigations have observed a positive correlation between gene body DNAhm and cis-gene expression. While some of these studies have observed an association between age and global DNAhm, none have investigated region-specific age-related DNAhm in human blood samples. In this study, we investigated DNAhm and gene expression in PBMCs of 10 young and 10 old, healthy female volunteers. Thousands of regions were differentially hydroxymethylated in the old vs. young individuals in gene bodies, exonic regions, enhancers, and promoters. Consistent with previous work, we observed directional consistency between age-related differences in DNAhm and gene expression. Further, age-related DNAhm and genes with high levels of DNAhm were enriched for immune system processes which may support a role of age-related DNAhm in immunosenescence.

The role of DNA methylation and hydroxymethylation in immunosenescence.

A healthy functioning immune system is critical to stave off infectious diseases, but as humans and other organisms age, their immune systems decline. As a result, diseases that were readily thwarted in early life pose nontrivial harm and can even be deadly in late life. Immunosenescence is defined as the general deterioration of the immune system with age, and it is characterized by functional changes in hematopoietic stem cells (HSCs) and specific blood cell types as well as changes in levels of numerous factors, particularly those involved in inflammation. Potential mechanisms underlying immunosenescence include epigenetic changes such as changes in DNA methylation (DNAm) and DNA hydroxymethylation (DNAhm) that occur with age. The purpose of this review is to describe what is currently known about the relationship between immunosenescence and the age-related changes to DNAm and DNAhm, and to discuss experimental approaches best suited to fill gaps in our understanding.

Non-linear patterns in age-related DNA methylation may reflect CD4+ T cell differentiation.

DNA methylation (DNAm) is an important epigenetic process involved in the regulation of gene expression. While many studies have identified thousands of loci associated with age, few have differentiated between linear and non-linear DNAm trends with age. Non-linear trends could indicate early- or late-life gene regulatory processes. Using data from the Illumina 450K array on 336 human peripheral blood samples, we identified 21 CpG sites that associated with age (P<1.03E-7) and exhibited changing rates of DNAm change with age (P<1.94E-6). For 2 of these CpG sites (cg07955995 and cg22285878), DNAm increased with age at an increasing rate, indicating that differential DNAm was greatest among elderly individuals. We observed significant replication for both CpG sites (P<5.0E-8) in a second set of peripheral blood samples. In 8 of 9 additional data sets comprising samples of monocytes, T cell subtypes, and brain tissue, we observed a pattern directionally consistent with DNAm increasing with age at an increasing rate, which was nominally significant in the 3 largest data sets (4.3E-15

Global epigenomic reconfiguration during mammalian brain development.

DNA methylation is implicated in mammalian brain development and plasticity underlying learning and memory. We report the genome-wide composition, patterning, cell specificity, and dynamics of DNA methylation at single-base resolution in human and mouse frontal cortex throughout their lifespan. Widespread methylome reconfiguration occurs during fetal to young adult development, coincident with synaptogenesis. During this period, highly conserved non-CG methylation (mCH) accumulates in neurons, but not glia, to become the dominant form of methylation in the human neuronal genome. Moreover, we found an mCH signature that identifies genes escaping X-chromosome inactivation. Last, whole-genome single-base resolution 5-hydroxymethylcytosine (hmC) maps revealed that hmC marks fetal brain cell genomes at putative regulatory regions that are CG-demethylated and activated in the adult brain and that CG demethylation at these hmC-poised loci depends on Tet2 activity.